Heterogeneous Stress Research Articles

Simulation Study on Diffusion of Chlorides in Concretes with Load-Induced Heterogeneous Stress

Simulation Study on Diffusion of Chlorides in Concretes with Load-Induced Heterogeneous Stress

Recent advances in characterizing the crustal stress field and future applications of stress data: Perspectives from North America

Abstract The stress field controls patterns of crustal deformation, including which faults are likeliest to cause earthquakes or transmit fluids. Since the 1950s, maps of maximum horizontal stress ( S Hmax ) orientations have advanced dramatically, and the style of faulting (relative principal stress magnitudes) has recently been mapped in some regions as well. This perspectives paper summarizes developments in characterizing stress orientations and (relative) magnitudes, including new seismic and borehole methods, as well as progress in identifying the causes of stress variations. Despite these advances, adding far more spatiotemporal detail would allow geoscientists to address many of today's key challenges regarding natural hazards, energy development, and geodynamics. In particular, it is critically important to characterize stress heterogeneity at multiple scales while also recognizing the coherent variability of the stress field. The second part of the paper considers how more detailed stress datasets could prove essential to addressing some of the grand questions in geoscience, including deciphering the poorly understood feedbacks between crustal dynamics and surface processes, improving earthquake and eruption forecasts, and determining the origins and shared properties of plate boundaries.

Predicting mechanical failure of polycrystalline dual-phase nickel-based alloys by numerical homogenization using a phase field damage model

AbstractBrazing of nickel-based alloys plays a major role in the assembly of turbine components, e.g., abradable sealing systems. In a brazed joint of nickel-based alloys a composition of brittle and ductile phases can be formed if the brazing conditions are not ideal. This heterogeneous microstructure is a crucial challenge for predicting the damage behavior of a brazed joint. The initiation and evolution of microdamage inside of the brittle phase of a virtual dual-phase microstructure representing the material in a brazed joint is studied by means of numerical simulations. A phase field approach for brittle damage is employed on the microscale. The simulation approach is capable of depicting phenomena of microcracking like kinking and branching due to heterogeneous stress and strain fields on the microscale. No information regarding the initiation sites and pathways of microcracks is needed a priori. The reliability of calculating the effective critical energy quantities as a microstructure-based criterion for macroscopic damage is assessed. The effective critical strain energy density and the effective critical energy release rate are evaluated for single-phase microstructures, and the approach is transferred to dual-phase microstructures. The local critical strain energy density turns out to be better suited as a model input parameter on the microscale as well as for a microstructure-based prediction of macroscopic damage compared to a model employing the energy release rate. Regarding the uncertainty of the model prediction, using the effective critical energy release rate leads to a standard deviation which is five times larger than the standard deviation in the predicted effective critical strain energy density.

A Source and Ground Motion Study of the Veracruz Earthquakes of 29 October 2009 (Mw5.7) and 4 August 2021 (Mw4.8): Evidence of Strong Azimuthal Variation of Attenuation

We study two moderate earthquakes that occurred offshore the State of Veracruz, in the southwestern Gulf of Mexico, on 29 October 2009 (Mw 5.7) and 4 August 2021 (Mw 4.8). The former was located near the town of Alvarado and latter near the city of Veracruz. The events were well recorded by accelerographs and seismographs at local and regional distances. W-phase regional centroid moment tensor inversion reveals that they had reverse-faulting mechanism, similar to several other earthquakes in the southwesternGulf of Mexico. Of the seven focal mechanisms now available along the southwestern margin, two are strike slip and the rest are of thrust type, suggesting a heterogeneous stress regime. We take advantage of local and regional recordings produced by these two earthquakes to study the characteristics of the ground motion. Source spectra computed at each station separately (withoutcorrecting for the site effect), assuming a reasonable geometrical spreading and Q = 141f 0.63, show remarkably high variability due to difference in path and local site effects. The geometric mean apparent source spectrum (source spectrum including site effects) of both earthquakes may be modeled by an ω2 -Brune source model with a stress drop, Δσ, of 40 MPa. These source spectra, along with the application of stochastic method, yield peak ground acceleration (PGA) and velocity (PGV) as a function of distance in general agreement with the observations. Of greater practical importance is the ground motion at sedimentary sites in the city of Veracruz and at the Laguna Verde Nuclear Power Plant (LVNPP) site, especially from a postulated Mw 6.5 earthquake which is a reasonable scenario event for the region. For the city of Veracruz and LVNPP we estimate site effect with respect to the ω2-Brune source with Δσ of 2 MPa. There is some support for this Δσ. We apply both stochastic and empirical Green’s functions (EGF) techniques in the estimation of the ground motion. The recording of the 2021 earthquake is taken as the EGF, and Δσ of the EGF and the target event are assumed to be the same and equal to 2 MPa. The predicted PGA and PGV at sedimentary sites in the city of Veracruz and LVNPP above the hypocenter (depth = 20 km) from the postulated Mw 6.5 earthquake are 0.2 g and 10 cm/s and 0.18 g and 3 cm/s, respectively. These results are preliminary as they are based on several assumptions.

Space Geodetic Insights to the Dramatic Stress Rotation Induced by the February 2023 Turkey‐Syria Earthquake Doublet

AbstractThe February 2023 Turkey‐Syria Earthquake doublet ruptured multiple segments of the East Anatolian Fault (EAF) Zone. Dominating seismicity focal mechanism shifted dramatically from strike‐slip to normal‐faulting after the doublet. To better understand this shift, here we derived a comprehensive 3D co‐seismic displacement field and performed the stress analysis. Abundant space geodetic data were used to generate high‐resolution 3D surface displacement, which provide tight constraints on fault geometry, slip distribution and stress field. Together with stress inversion from aftershock focal mechanisms, we show that the principal stress direction rotation in the region with the most normal‐faulting aftershocks is the staggering 29°. The induced heterogenous stress may explain the shift of the dominant focal mechanism toward normal faulting. We suggest that the extensional horsetail splay faults, likely formed through geologic time scale related to the releasing bend on the EAF, are the hosts of most of the normal faulting aftershocks.

Effect of pore on the deformation behaviors of NiTi shape memory alloys: A crystal-plasticity-based phase field modeling

Porous NiTi shape memory alloys (SMA) have been developed into important biomedical and damping materials. However, the effect of pore on the martensitic microstructure evolution, functional properties, and plastic deformation has not been fully clear yet. In this work, to reveal the functional properties including superelasticity (SE), elastocaloric effect (eCE), one-way shape memory effect (OWSME), and stress-assisted two-way shape memory effect (SATWSME), as well as plastic deformation of porous NiTi SMAs, the polycrystalline SMA system is considered as a composite containing the grain-interior (GI) and grain-boundary (GB) phases, a thermo-mechanically coupled, grain-size-dependent and crystal-plasticity-based phase field model is proposed for the GI phase, and an elasto-viscoplastic constitutive model is established for the GB phase. The effects of porosity, pore diameter, pore aspect ratio, pore orientation, and pore distribution, as well as the interaction between pore and grain size, are comprehensively revealed and discussed. The simulations demonstrate that the pores cause a heterogeneous stress field, and the interaction among pores and that between pores and grain boundaries can enhance such heterogeneity, which can promote the initiation of martensitic transformation and reorientation. Moreover, the geometrical constraints caused by the pores can significantly affect the morphology and evolution of martensite microstructures. However, due to the diverse microstructure evolutions during SE, OWSME, and SATWSME, their interactions with the pores are quite different, resulting in various dependencies of each functional property and plastic deformation on the pore, and such dependencies are strongly influenced by porosity, pore diameter, pore aspect ratio, pore orientation, and pore distribution. The simulated results in this work are helpful to comprehensively and deeply understand the pore effect and its microscopic mechanism.

On the Unloading‐Induced Fault Reactivation: The Effect of Stress Path on Failure Criterion and Rupture Dynamics

AbstractFault reactivations induced by deep excavation can pose significant challenges to underground construction or resource extraction. Laboratory experiments on rock faults demonstrate that unloading‐induced fault reactivations obey the Coulomb failure criterion derived from loading‐induced events. However, the effect of stress path during unloading on the failure criterion and rupture dynamics of fault reactivations remains poorly understood. Here, we present findings from a series of laboratory experiments aimed at elucidating the effect of the unloading path on the failure criterion and rupture dynamics of fault reactivations. We conducted experiments under various stress conditions, examining two cases of unloading paths. In Case I, we unloaded the minimum principal stress, while in Case II, the maximum principal stress was unloaded. Strain gauges and high‐speed photography were employed to capture the transient dynamic rupture process. Our investigations have yielded new insights into the effect of unloading path on the rupture dynamics when the fault is reactivated. In Case I, we observed fault reactivations resembling those loading‐induced events characterized by forward sliding. Conversely, in Case II, fault reactivations associated with stress reversal produce mild reversed sliding with lower stress drop and rupture velocity. Furthermore, we find that there is a remarkable reduction in static friction for reversed sliding, indicating that the failure criterion for fault reactivation is influenced by the stress path. We demonstrate that enhanced stress heterogeneity, caused by stress reversal, serves as a mechanism for reduced static friction. These findings contribute to our understanding of the mechanisms underlying fault reactivations, particularly those involving reversed sliding.

The combined effects of a heterogeneous porosity distribution and stress gradient on the high cycle fatigue behavior of high pressure die cast AlSi9Cu3

In order to expand the number of possible applications of High Pressure Die Cast (HPDC) components using the AlSi9Cu3 alloy, the automotive industry is giving greater attention to the fatigue design of components produce by this process-material combination. The process is well known for introducing a heterogeneous distribution of casting defects. This study explores the impact of the heterogeneous porosity distribution and stress gradients on the high cycle fatigue behavior of HPDC AlSi9Cu3 components. Fatigue tests under fully-reversed tension–compression and fully-reversed bending loads are used to investigate the influence of porosity and stress gradients on the fatigue behavior, by machining different specimen thicknesses. The results demonstrate that stress gradients and the heterogeneous porosity distribution determine the material’s fatigue response, including both the fatigue strength and crack initiation mechanisms. The Murakami parameter and local stress level at the critical defect make it possible to interpret the gradient effects. The nature of the entrapped gas in internal defects may influence the in-bulk crack propagation behavior. The comparison of fracture mechanics models (El Haddad, Murakami, Kitagawa-Takahashi) make it possible to highlight the link between the El Haddad’s parameter to parameters that are more microstructurally based.

Effects of Different Water Stress Levels, Heterogeneity and Location on Berry Phytochemical Properties in an Organic and Conventional Vineyard (Vitis vinifera cv. Cabernet-Sauvignon)

The research was conducted in two vineyards, one certified as organic and the other following conventional practices. The experiment, designed with Split-Plot Experimental Design based on stress levels, was divided into two main plots, Organic and Conventional, and each of these plots was further divided into two subplots. The predawn leaf water potential results categorized the vines into two groups: those with values lower than -8 MPa and those greater than -8 MPa, which were labeled as Dryland-shallow soil and Baseland-deep soil, respectively, based on the location and soil type. During the harvest, grape clusters were collected and classified into three groups based on berry diameter (10mm-12mm, 12mm-14mm, 14mm-16mm). The results indicated that the 10mm-12mm berry size group generally exhibited the desired characteristics across all evaluated criteria. The total anthocyanin and total tannin content were found to be higher in the 10mm-12mm berries from vines experiencing moderate stress (Stress 1), regardless of their location. Additionally, the Dryland-shallow soil condition showed higher tannin content. On the other hand, grapes from high-stress vines displayed lower antioxidant values. The total polyphenol index content was higher in the organic vineyard. Based on the findings, it was suggested that to obtain high phytochemical compounds from Cabernet-Sauvignon grapes in the Tekirdağ region, cultivation should be carried out under Dryland-shallow soil conditions, where the predawn leaf water potential can drop as low as -0.8 MPa during the period between veraison and harvest. Moreover, using berries with size between 10mm-12mm might be suitable for this purpose.

Optimizing extreme limited entry design for lateral heterogeneous reservoirs

Extreme limited entry (XLE) method can enhance hydraulic fracture uniformity by elevating wellbore pressure to overcome reservoir heterogeneity. However, it increases surface pressure and decreases hydraulic fracturing efficiency. In this paper, the effects of stress shadow and heterogeneity of the lateral breakdown pressure on fracture propagation are numerically studied, from which an optimization method is proposed for the XLE design with the minimal requirement of fracturing pressure. Initially, the unconventional fracture model is applied to analyze fracture propagation in formations with homogeneous geomechanical parameters and various perforation designs. Fracture length distribution is correlated with perforation friction, revealing insights into the influence of geomechanical parameters and fracture spacing. Then, the uniformity of fracture propagation is assessed for different XLE designs by varying lateral breakdown pressures. Simulation results are utilized to develop a cross-plot of fracture length distributions vs normalized perforation frictions for optimizing an XLE design. In homogeneous reservoirs, decreasing perforation spacing and increasing perforation hole erosion rate slow down the reduction in the fracture length distribution span with perforation friction, necessitating a higher perforation friction for uniform fractures. In all simulated cases, a similar trend is observed when the distribution of fracture lengths is plotted vs the normalized perforation friction to the difference of breakdown pressure within one fracturing stage. It is recommended that this normalized value be maintained above 2 in an XLE design. This new plot makes the conventional horn-shaped plot more applicable to reservoirs with heterogeneous lateral breakdown pressures.

Dendrite Growth—Microstructure—Stress—Interrelations in Garnet Solid‐State Electrolyte

AbstractThis study illustrates how the microstructure of garnet solid‐state electrolytes (SSE) affects the stress‐state and dendrite growth. Tantalum‐doped lithium lanthanum zirconium oxide (LLZTO, Li6.4La3Zr1.4Ta0.6O12) is synthesized by powder processing and sintering (AS), or with the incorporation of intermediate‐stage high‐energy milling (M). The M compact displays higher density (91.5% vs 82.5% of theoretical), and per quantitative stereology, lower average grain size (5.4 ± 2.6 vs 21.3 ± 11.1 µm) and lower AFM‐derived RMS surface roughness contacting the Li metal (45 vs 161 nm). These differences enable symmetric M cells to electrochemically cycle at constant capacity (0.1 mAh cm−2) with enhanced critical current density (CCD) of 1.4 versus 0.3 mA cm−2. It is demonstrated that LLZTO grain size distribution and internal porosity critically affect electrical short‐circuit failure, indicating the importance of electronic properties. Lithium dendrites propagate intergranularly through regions where LLZTO grains are smaller than the bulk average (7.4 ± 3.8 µm for AS in a symmetric cell, 3.1 ± 1.4 µm for M in a half‐cell). Metal also accumulates in the otherwise empty pores of the sintered compact present along the dendrite path. Mechanistic modeling indicates that reaction and stress heterogeneities are interrelated, leading to current focusing and preferential plating at grain boundaries.

Can Geometrical Barrier Explain the Mw 7.8 Earthquake in Southern Türkiye on February 2023?

Abstract Two large-magnitude earthquakes hit southern Türkiye on February 2023. The first, Mw 7.8 strike-slip earthquake generated a rupture of 300 km section along the ∼600 km long East Anatolian fault (EAF). Here, we present an analytical solution using perturbation theory for the static stress field near the EAF induced by the fault geometry and the tectonic loading before these earthquakes. By applying the Coulomb failure criterion, we show that a large stress barrier is developed around the segment that ruptured in the first earthquake. Considering stress field conditions that are associated with left-lateral strike-slip on the fault, we demonstrate how the barrier location is mostly determined by the fault geometry, while its magnitude is sensitive to the background stress value and direction. We further show that the elastic energy around the fault increases to maximum values near the barrier region and decreases away from it. Therefore, we suggest that the high magnitude and the associated long rupture of the earthquake were strongly influenced by the static stress heterogeneity generated by the fault geometry.

Fault roughness controls injection-induced seismicity

Surface roughness ubiquitously prevails in natural faults across various length scales. Despite extensive studies highlighting the important role of fault geometry in the dynamics of tectonic earthquakes, whether and how fault roughness affects fluid-induced seismicity remains elusive. Here, we investigate the effects of fault geometry and stress heterogeneity on fluid-induced fault slip and associated seismicity characteristics using laboratory experiments and numerical modeling. We perform fluid injection experiments on quartz-rich sandstone samples containing either a smooth or a rough fault. We find that geometrical roughness slows down injection-induced fault slip and reduces macroscopic slip velocities and fault slip-weakening rates. Stress heterogeneity and roughness control hypocenter distribution, frequency-magnitude characteristics, and source mechanisms of injection-induced acoustic emissions (AEs) (analogous to natural seismicity). In contrast to smooth faults where injection-induced AEs are uniformly distributed, slip on rough faults produces spatially localized AEs with pronounced non-double-couple source mechanisms. We demonstrate that these clustered AEs occur around highly stressed asperities where induced local slip rates are higher, accompanied by lower Gutenberg-Richter b-values. Our findings suggest that real-time monitoring of induced microseismicity during fluid injection may allow identifying progressive localization of seismic activity and improve forecasting of runaway events.

Exploring the role of Type-II residual stresses in a laser powder bed fusion nickel-based superalloy using measurement and modeling

Far-field high-energy X-ray diffraction microscopy (ff-HEDM) and the crystal plasticity finite element method (CPFEM) are used to investigate the role of grain-scale (Type-II) residual stresses on the fatigue life of additively manufactured (AM) Inconel alloy 625 (IN-625). Grain-averaged orientations, centroids, and residual elastic strain tensors from ff-HEDM data are used to instantiate a crystal plasticity model to simulate the effect of residual stresses at the grain scale. Simulation results indicate that the presence of tensile residual strains increase stress localization and heterogeneity within grains, triggering an earlier onset of plasticity. A microscale fatigue indicator parameter (FIP) is computed to model the impact of these residual strains on the cycles to fatigue crack nucleation. The crack nucleation model, based on the computed FIPs, predicts a significant reduction in the number of cycles for fatigue crack nucleation for mid- and high-cycle fatigue due to the residual strain induced localization, while the residual strains have minimal impact on low-cycle fatigue life.

Neopolyploidy increases stress tolerance and reduces fitness plasticity across multiple urban pollutants: support for the “general-purpose” genotype hypothesis

Abstract Whole-genome duplication is a common macromutation with extensive impacts on gene expression, cellular function, and whole-organism phenotype. As a result, it has been proposed that polyploids have “general-purpose” genotypes that perform better than their diploid progenitors under stressful conditions. Here, we test this hypothesis in the context of stresses presented by anthropogenic pollutants. Specifically, we tested how multiple neotetraploid genetic lineages of the mostly asexually reproducing greater duckweed (Spirodela polyrhiza) perform across a favorable control environment and 5 urban pollutants (iron, salt, manganese, copper, and aluminum). By quantifying the population growth rate of asexually reproducing duckweed over multiple generations, we found that across most pollutants, but not all, polyploidy decreased the growth rate of actively growing propagules but increased that of dormant ones. Yet, when considering total propagule production, polyploidy increased tolerance to most pollutants, and polyploids maintained population-level fitness across pollutants better than diploids. Furthermore, broad-sense genetic correlations in growth rate among pollutants were all positive in neopolyploids but not so for diploids. Our results provide a rare test and support for the hypothesis that polyploids are more tolerant of stressful conditions and can maintain fitness better than diploids across heterogeneous stresses. These results may help predict that polyploids may be likely to persist in stressful environments, such as those caused by urbanization and other human activities.

Agent-based modeling of stress anisotropy driven nematic ordering in growing biofilms.

Living active collectives have evolved with remarkable self-patterning capabilities to adapt to the physical and biological constraints crucial for their growth and survival. However, the intricate process by which complex multicellular patterns emerge from a single founder cell remains elusive. In this study, we utilize an agent-based model, validated through single-cell microscopy imaging, to track the three-dimensional (3D) morphodynamics of cells within growing bacterial biofilms encased by agarose gels. The confined growth conditions give rise to a spatiotemporally heterogeneous stress landscape within the biofilm. In the core of the biofilm, where high hydrostatic and low shear stresses prevail, cell packing appears disordered. In contrast, near the gel-cell interface, a state of high shear stress and low hydrostatic stress emerges, driving nematic ordering, albeit with a time delay inherent to shear stress relaxation. Strikingly, we observe a robust spatiotemporal correlation between stress anisotropy and nematic ordering within these confined biofilms. This correlation suggests a mechanism whereby stress anisotropy plays a pivotal role in governing the spatial organization of cells. The reciprocity between stress anisotropy and cell ordering in confined biofilms opens new avenues for innovative 3D mechanically guided patterning techniques for living active collectives, which hold significant promise for a wide array of environmental and biomedical applications.

Неотектоника и геодинамика зоны сочленения Северо-Западного Кавказа и Черноморской впадины (район г. Геленджика)

Relevance of the work. In the junction zone of the neotectonic structures of the North-West Caucasus and the Black Sea depression, structures of different types in strike and morphological-kinematic type are developed. Determining the structural-geomorphological and neotectonic zoning, including modern ones, and its connection with the sources of tectonic forces and stresses is relevant. Such studies are necessary to ensure the geodynamic safety of engineering structures in the city of Gelendzhik, the territory of which is subject to negative geological processes. The purpose of the research is to identify the structural-geomorphological, neotectonic modern zoning with an assessment of the modern geodynamic conditions of its formation in connection with ensuring the safety of engineering drainage routes in the Gelendzhik region. The research method is structural-geomorphological, with the construction of profiles and geomorphological steps with an assessment of their age and deformations. The work used neotectonic and structural-geodynamic approaches aimed at identifying structures and identifying the sources of stress and deformation fields. The research methodology is the recognition of heterogeneity, different ranks and different ages of stresses and deformations. As a result of the work, erosion-denudation stages of different ages and marine terrace levels and neotectonic structures were identified with an assessment of modern geodynamic conditions. For the first time, a cross-structural plan has been established, in which all-Ural and anti-Caucasian structures intersect. Their origin is associated with the superposition of heterogeneous forces and stresses. Compression stress is generated in connection with collision processes of lithospheric plates and microplates (blocks), tension is generated in connection with the activity of the deep East Black Sea geodynamic system. It considers as a mechanism the upper mantle convection directed from the Black Sea depression towards the orogen of the North-Eastern Caucasus. The superposition of faults of different ages and different types on top of each other causes increased activity in them. Late orogenic faults, consistent with early Alpine (early orogenic) reverse faults, are characterized by increased displacement amplitudes and wide dislocation zones. It is dangerous for engineering routes to cross such gaps. The pipeline at the intersection with the zone of active ruptures will be equipped with special hinge-type devices to prevent dangerous deformations in it.

NUMERICAL INVESTIGATION OF THE VALVE MODEL FOR PRIMARY SUTURELESS IMPLANTATION

HighlightsThe presented numerical model of the bioprosthesis support frame for sutureless fixation needs to be optimized to reduce the amplitude of the von Mises stresses.The presented setup of numerical modeling has been validated using the study of commercial analogs. AbstractAim. To assess the stress-strain state of three computer models of TAVI prosthesis support frames during crimping and shaping to the delivery configuration.Methods. The study included three models of stent-like support frames for balloon-expandable devices, which are intended to serve as a foundation for the subsequent development of a domestic aortic valve prosthesis: two models of commercial bioprostheses and one experimental. The objects were evaluated numerically under the conditions of stress-strain state that arises in the stent support frames, simulating two loads: crimping and shaping to the delivery configuration. The study was conducted using the numerical modeling complex Abaqus/CAE (“Dassault Systemes”, France). The key indicators for evaluation were the von Mises stress, as a strength criterion, and its distribution over the frame; the presence and proportion of elastic recoil.Results. It was shown that all samples are capable of reaching the required diameter without excessive material. A quantitative investigation of the Von Mises stress showed that commercial models display amplitudes below the material`s strength limit (892.4 and 916.8 MPa), whereas the proprietary model exceeds this limit, reaching 991.4 MPa, requiring geometry optimization. Shaping to the delivery configuration indicates that all models provide safe expansion up to 26 mm, with a Von Mises stress level in the range of 882.4–914.1 MPa, which is below the strength limit of the cobalt-chrome alloy. Moreover, we have noted heterogeneous stress distribution, with concentration in the lamella junction areas.Conclusion. Thus, it has been demonstrated that numerical modeling and the finite element method can be effectively applied to assess the stress-strain state of sutureless prostheses. Geometry optimization and further development of this project – the development of a Russian minimally invasive aortic valve prosthetic system, may contribute towards increasing the accessibility of this treatment method.

Integration of geological and well log and geomechanical modeling toward a safe mud weight for wellbore stability: A case study of Ruby field, Blocks 01&02, Cuu Long basin

There is a high density of drilled wells in Blocks 01&02 which is a part of the Cuu Long basin. The geological settings of this area have undergone a complicated evolution, resulting in heterogeneous lithology and variable stresses. Consequently, drilling activities in the Ruby oil field have faced numerous challenges, such as stuck pipe, gas kick, mud loss, and even lost well, etc. leading to substantial time and financial setbacks. It is essential to thoroughly understand the geomechanical characteristics of Ruby field to ensure safe operation and optimize drilling. In this paper, several downhole geophysical data sets such as wireline logging data as gamma ray (GR), density log (RHOB), neutron (NPHI), compression/shear sonic travel time (DTC/DTS), image formation log (FMI), pressure test (MDT) and leak of test (LOT, XLOT) are used to construct the 1D geomechanical model. The objective of this modelling is to compute parameters of pore pressure, vertical stress, horizontal stress, elastics properties, rock strengths. Then, these parameters are used to analyze the wellbore stability and to recommend appropriate drilling mud weights for the wells under study. This analysis can subsequently be extended to cover the entire Ruby oil field for the future drilling operations to enhance overall safety and efficiency.

The role of interface geometry and appendages on the mesoscale mechanics of the skin.

The skin is the largest organ in the human body and serves various functions, including mechanical protection and mechanosensation. Yet, even though skin's biomechanics are attributed to two main layers-epidermis and dermis-computational models have often treated this tissue as a thin homogeneous material or, when considering multiple layers, have ignored the most prominent heterogeneities of skin seen at the mesoscale. Here, we create finite element models of representative volume elements (RVEs) of skin, including the three-dimensional variation of the interface between the epidermis and dermis as well as considering the presence of hair follicles. The sinusoidal interface, which approximates the anatomical features known as Rete ridges, does not affect the homogenized mechanical response of the RVE but contributes to stress concentration, particularly at the valleys of the Rete ridges. The stress profile is three-dimensional due to the skin's anisotropy, leading to high-stress bands connecting the valleys of the Rete ridges through one type of saddle point. The peaks of the Rete ridges and the other class of saddle points of the sinusoidal surface form a second set of low-stress bands under equi-biaxial loading. Another prominent feature of the heterogeneous stress pattern is a switch in the stress jump across the interface, which becomes lower with respect to the flat interface at increasing deformations. These features are seen in both tension and shear loading. The RVE with the hair follicle showed strains concentrating at the epidermis adjacent to the hair follicle, the epithelial tissue surrounding the hair right below the epidermis, and the bulb or base region of the hair follicle. The regions of strain concentration near the hair follicle in equi-biaxial and shear loading align with the presence of distinct mechanoreceptors in the skin, except for the bulb or base region. This study highlights the importance of skin heterogeneities, particularly its potential mechanophysiological role in the sense of touch and the prevention of skin delamination.